Path Entanglement of Continuous- Variable Quantum Microwaves 
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Path entanglement constitutes an essential resource in quantum information and communication 
protocols. Here, we demonstrate frequency-degenerate entanglement between continuous- variable 
quantum microwaves propagating along two spatially separated paths. We combine a squeezed and 
a vacuum state using a microwave beam splitter. Via correlation measurements, we detect and 
quantify the path entanglement contained in the beam splitter output state. Our experiments open 
the avenue to quantum teleportation, quantum communication, or quantum radar with continuous 
variables at microwave frequencies. 



PACS numbers: 03.67.Bg, 03.65.Ud, 42.50.Dv, 85.25.-j 

Fascinatingly, quantum mechanics allows for a com- 
pound system to have a common description while, at 
the same time, no individual states can be ascribed 
to its subsystems^. The presence of entanglement be- 
tween spatially separated systems is a necessary condi- 
tion for what Einstein called "spooky action at a dis- 
tance"^!: the contradiction between quantum mechan- 
ics and local realisrrP^. Furthermore, entanglement is 
at the heart of quantum communication and informa- 
tion processing technologies, which promise significant 
performance gains over classical protocol d 1 * 4 * 5 *. Conse- 
quently, entanglement has been extensively explored in 
atomic physics and quantum optics^™. In these in- 
vestigations, optical frequencies were preferred over mi- 
crowaves because the higher photon energies facilitate 
practical applications. However, since the late 1990s, mi- 
crowave technology has evolved rapidly in both industry 
and science. For one thing, classical microwave fields 
have become an indispensable tool in mobile communi- 
cation. For another, a promising direction towards scal- 
able quantum information processing has appeared with 
the advent of superconducting microwave quantum cir- 
cuit*^. Despite some decoherence issues, these systems 
provide unprecedented light-matter coupling strengths 
due to their large effec tive dipole moments and field en- 
hancement effects^HU, As a consequence, standing- wave 
fields in transmission line resonators were shown to act 
as a sho rt-ran ge quantum bus between superconduct- 
ing qubits^^ and various gates were implemented^^. 
For microwave quantum communication, however, prop- 
agating fields are required. As a first step in this di- 
rection, early experiments demonstrated tomography of 
weak thermal states^, coherent states^, and single pho- 



tons^. Next, continuous- variable states generated by 
Josephson parametric devices were reconstructed^. Very 
recently, such devices have permitted to investigate two- 
mode squeezin^M^. An important aspect of these ex- 
periments is the understanding they provide regarding 
entanglement. In order to be a resource in quantum com- 
munication protocols, it must occur between spatially 
separated subsystems^. Furthermore, a strict proof of en- 
tanglement requires the entangler and the detector to be 
based on independent experimental techniques. In this 
work, we make a significant step beyond previous efforts 
and demonstrate path entanglement in the microwave 
regime, respecting both criteria mentioned above. Our 
experiments follow the spirit of the quantum-optical re- 
alization^ of the original Einstein- Podolsky- Rosen (EPR) 
paradox^. As shown in Fig. [I] we combine a vacuum and 
a squeezed vacuum state in a hybrid ring microwave beam 
splitter 2 - 4 ^ acting as an entangling device. Its two output 
ports hold a continuous- variable state which is frequency- 
degenerate and entangled with respect to the two prop- 
agation paths. Along these paths, the entanglement can 
be conveniently distributed to two parties requiring it 
for any suitable quantum communication protocol. In 
our experiments, we first reconstruct the squeezed in- 
put state by means of dual-path tomograph}^, which 
assumes knowledge of the beam splitter relations. Next, 
we reconstruct the moments of the output state after the 
beam splitter by treating the latter as a black box and 
calibrating against a known stated In this reference- 
state method (see Supplementary), we only assume that 
independent vacuum states are produced in each out- 
put path when vacuum is incident at both input ports. 
From the moments reconstructed in this way, we build 
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FIG. 1. Layout of the experiment. The microwave beam splitter acts as entangling device (green eight-shaped structure). The 
blue-and- red arrows denote the path-entangled state. In the cross correlation detector, the oscilloscope symbols denote the 
noisy amplification, down- conversion, and digitizing of the in-phase (ii.,2) and quadrature (Qi,2) components of the output 
signals. The computer symbolizes the final numerical data processing, partly done by an FPGA logic. 



a witness matrix which proves the existence of path en- 
tanglement independently of the detailed nature of our 
output stated Since in practice the data shows that 
our states are Gaussian, we finally quantify the degree 
of entanglement by means of the negativity^- 7 . The re- 
sult of this analysis agrees with what we expect for our 
squeezed input state. We note here that for bipartite 
single- mode Gaussian entanglement, as it is relevant in 
our case, entanglement implies nonlocality 1 5 . All in all, 
our results show that we have realized the main building 
block for microwave quantum teleportation and commu- 
nication protocols. 

The generation of the input states for the beam splitter 
is straightforward. The vacuum is realized with a com- 
mercial 50 (^-termination at 40 mK acting as a broadband 
black-body emittei^. The squeezed state is produced us- 
ing a particular superconducting circuit, the flux-driven 
Josephson parametric amplifier (JPA) 28 . In this device 
two Josephson junctions form a nonlinearity which can 
be modulated ("pumped") at gigahertz frequencies to 
achieve a parametric effect. The JPA box is stabilized to 
50 mK. A thermal state emitted by an attenuator, whose 
temperature can be varied from 50 — 800 mK, can be fed 
into the JPA. Our cross correlation detector is based on 
the insight that for microwave signals off-the-shelf high- 
gain low-noise linear amplifiers are available rather than 
efficient single photon counters. We connect one amplifi- 
cation path to each output port of the beam splitter. At 
room temperature we record the in-phase and quadra- 
ture components, and of the amplified signals. 

The averaged moments (I1I2QTQ2) are computed for 
j + /c + m + n<4 and j, fc, m, n G No in real time using a 
field programmable gate array (FPGA) logic. Further 
details can be found in the Supplementary. 



As a first test of our setup, we perform dual-path recon- 
structions of the Wigner function for known input states. 
Here, we exploit the fact that the noise contributions of 
the two amplification paths are independent, while the 
split signals are correlatecP^(see Supplementary). We re- 
construct vacuum fluctuations and coherent states (dis- 
placed vacuum), both at a frequency fo = 5.637 GHz. Be- 
cause of narrow-band filtering, we approximate the vac- 
uum and thermal states as single-mode fields. The results 
shown in Fig. [2^a) exhibit a very good phase control for 
the coherent state. In addition, we find a small thermal 
contribution of 0.097±0.007 photons above the vacuum 
level which can be due to a small thermal population or 
other experimental imperfections. In the next step, we 
generate a squeezed state by pumping the JPA. For a sig- 
nal gain of 10 dB and a phase of 45°, the reconstructed 
Wigner function is shown in Pig. j2^b). An analysis of the 
reconstructed signal moments reveals that, at the input 
of the beam splitter, the state generated by the JPA is 
squeezed by 4.9±0.2 dB below the vacuum level and con- 
tains 8. 72 ±0.05 photons. Furthermore, the product of 
the standard deviation of the squeezed quadrature with 
that of its orthogonal, enlarged one, is 3.45±0.07 times 
larger than the variance of the ideal vacuum. In other 
words, we can model the state as one created by an 
ideal squeezer acting on an effective thermal field with 
1.22±0.04 photons. This thermal field contains the com- 
bined effects of losses and the small thermal population 
found in the experimental vacuum. Again, we notice 
good control of the phase. It is noteworthy to mention 
that the amount of squeezing quoted above is mainly lim- 
ited by cable losses and not by the JPA itself. 

After characterizing the input fields of the beam split- 
ter, we now turn to its outputs. With the reference-state 
method, we build an entanglement witness matrix from 
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FIG. 2. Dual-path reconstruction of various states incident at the "squeezed state input port" of the beam splitter, p and q are 
dimensionless variables spanning the phase space, (a) JPA pump off. Reconstruction of the vacuum and of displaced vacuum 
states (coherent states, 8.80±0.01 photons, eight different phase values). All nine Wigner functions are superposed, (b) JPA 
pump on. Squeezed state for 10 dB JPA signal gain at 45°. Inset: 1/e contours of the ideal vacuum (blue), the experimental 
vacuum (green) displayed in panel (a), and the squeezed state (red). 



the reconstructed moments. Our witness reliably dis- 
tinguishes between "separable outputs" for the vacuum 
state and "path entangled outputs" for the squeezed state 
input. Next, we analyze the third and fourth order cu- 
mulants and find them to be small for JPA signal gains 
up to 10 dB. Since this is a strong indication for Gaus- 
sian states, we explore the path entanglement generated 
in our setup quantitatively via the negativity A/* ou t- For 
positive values, Af ou t describes the degree of entangle- 
ment produced between the beam splitter output paths 
(see Supplementary). In the limit of low JPA signal gain, 
Fig.[3ja) shows how A/" ou t becomes suppressed when send- 
ing more and more thermal photons into the JPA. At 
some point, the JPA cannot squeeze the incoming field 
below the vacuum anymore and the output state is no 
longer entangled. For constant temperature, Fig. [3^b) 
shows how A/"out increases with increasing signal gain from 
zero to a value A/" ou t,max = 0.55±0.04 at 10 dB signal gain. 
This behavior is in good agreement with the negativity 
A/*caic calculated from the dual-path reconstructed input 
state. Again, we observe a suppression for large thermal 
fields sent into the JPA. Our results confirm the expecta- 
tion 29 that the degree of squeezing at the beam splitter 
input determines the amount of entanglement generated 
between the output paths. However, since A/caic is gener- 
ally slightly lower than A/" ou t 5 we conclude that either the 
dual-path reconstruction underestimates the squeezing at 
the beam splitter input or the reference-state method ig- 
nores a small amount of spurious classical correlations be- 
tween the two paths. Both effects are consistent with the 
data shown in Fig. |3Fa), where at constant signal gain, 



the curve measured with the reference-state method at 
the beam splitter output converges for high temperatures 
to that calculated from the dual-path reconstructed in- 
put state. We finally note that the path-entangled state 
is expected to be a two-mode squeezed state with two ad- 
ditional local squeezing operations applied to itP^l Since 
local operations do not change the amount of entangle- 
ment, the negativity A/" ou t,max = 0.55±0.04 implies that 
the two-mode squeezed state before the two local oper- 
ations would have a variance squeezed by 3.2 dB below 
that of the two-mode vacuum. 

In summary, we present clear evidence for path 
entanglement generated by combining two frequency- 
degenerate continuous- variable microwave fields, the vac- 
uum and the squeezed vacuum, in a beam splitter. For 
an input state squeezed 4.9±0.2dB below the vacuum, 
we observe a maximum negativity A/" ou t,max = 0.55±0.04 
at 10 dB JPA signal gain. Our experiments bring the 
exciting quantum physics of entangled propagating elec- 
tromagnetic fields to the technologically highly attrac- 
tive microwave domain. In this way, they open up new 
and exciting perspectives towards microwave quantum 
teleportation, quantum communication, and quantum 
radaM 
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FIG. 3. Quantitative analysis of the path entanglement generated in our experiments. The negativities A/" ou t, A/" ca ic are the 
maxima of the corresponding negativity kernels Afout, A/'caic and 0. Circular symbols: A/" ou t data at the beam splitter output. 
Square symbols: A/'caic calculated from the reconstructed input state. The lines are guides to the eye. (a) Negativity kernel 
versus attenuator temperature (color code) at 1 dB signal gain. For the data points in the shaded area, the witness matrix 26 
confirms entanglement, (b) Negativity kernel versus the JPA signal gain. The blue (red) curves are recorded at 50 mK (573 mK). 
Grey point: negativity of the reference state, assumed to be zero. 
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Supplementary 



I. EXPERIMENTAL DETAILS 

A. The flux driven Josephson parametric amplifier 

We generate squeezed states using a flux driven 
Josephson parametric amplifier^ (JPA). Micrographs 
of the device used in our experiments are shown 
in Figs. ^TJa)-(d), a circuit diagram is displayed in 
Fig. £[lje). In order to achieve a parametric effect, the 
resonance frequency f^ c of a quarter wavelength super- 
conducting coplanar waveguide resonator (external qual- 
ity Qext = 312) is modulated at 2/dc- To this end, a dc 
superconducting quantum interference device (SQUID) 
- a superconducting loop interrupted by two Josephson 
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Fig. SI. Flux driven Josephson parametric amplifier used 
in our experiments, (a) Sample chip, (b) Zoom-in of the 
region marked with the red rectangle in panel (a), showing the 
coupling capacitor C c - (c) Pump line and dc SQUID. Zoom- 
in of the region marked with the blue rectangle in panel (a), 
showing the pump line and the dc SQUID, (d) Zoom-in of the 
region marked with the green rectangle in panel (c), showing 
the dc SQUID, (e) Circuit diagram. 



junctions - is inserted between the center conductor and 
the ground plane at the shorted end of the resonator. 
Because the SQUID acts as a flux-tunable inductor, the 
resonance frequency of the resonator can be changed by 
applying an external magnetic field. Via an external coil, 
we first set a quasistatic bias corresponding to a JPA op- 
erating frequency of fo = 5.637 GHz (see Fig. ^2|. The 
fast modulation at 2/o, also referred to as th e pum p tone, 
is then applied in a pulsed fashion (see Sec. II A[ ) via an 
on-chip antenna. We determine the isolation between 
antenna and resonator to be at least 28 dB at the signal 
frequency fo. 
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Fig. S2. JPA resonance frequency /dc as a function of the 
applied dc flux $dc- Black symbols: data. Red line: fit. Blue 
dot: operating point /o = 5.637 GHz. 

Resonator and antenna are made of a 50 nm thick Nb 
film. At the contacts, 95 nm of gold on a 5nm tita- 
nium bonding layer are deposited on top. As substrate, 
we use thermally oxidized (300 nm) silicon with a thick- 
ness of 300 |Ltm. The dc SQUID is fabricated in the 
last step using aluminum technology and shadow evap- 
oration. The Al electrodes have a thickness of 50 nm 
each. From Fig. ^2j we estimate a Josephson cou- 
pling energy Ej jh = 650 GHz for each junction, where 
h = 6.63 x 10 -34 Js is the Planck constant. The sample 
chip is placed between two small alumina printed circuit 
boards inside a gold-plated copper box. 



B. The 180° hybrid ring microwave beam splitter 

Microwave beam splitters divide an input signal while 
possibly adding phases to their outputs. However, in 
order to be lossless, matched, and reciprocal, these de- 
vices must necessarily have four ports^. In a quantum 
mechanical picture, this implies that a second, possibly 
hidden input port is always presents 3 . A sketch of the 
180° hybrid ring microwave beam splitter used in this 
work is shown in Fig. ^3j It is a commercially fabri- 
cated device based on gold microstrip transmission lines 
on a dielectric substrate encased in a copper housing. Its 
functionality is best understood as follows: at a wave- 
length A, the signals incident at the input ports form 
an interference pattern in the ring with antinodes at the 
output ports and nodes at the input ports 32 . This was 
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Fig. S3. Schematic sketch of the 180° hybrid ring microwave 
beam splitter used in our experiments. The split signals from 
the JPA acquire a 180° phase shift with respect to each other, 
while the split signals from the 50 Q termination remain in 
phase. 



experimentally demonstrated in Ref. |3H where also typ- 
ical transmission characteristics are shown. The signals 
from the two input ports are evenly split and superposed 
in the two output ports. In this superposition, the com- 
ponent from one of the inputs acquires a 180° phase shift 
between the output ports, while that from the other one 
remains in phase. The center frequency of our hybrid 
ring is 5.75 GHz. Within our measurement bandwidth 
of 2 x B Wfiiter = 2 x 489 kHz centered at the JPA operat- 
ing frequency /o (see also Fig. ^8|, our hybrid ring still 
has a coupling of 3.5 dB between input and output ports 
and an isolation of at least 38 dB between any two input 
or output ports. The magnitude imbalance between the 
two output ports is only 0.03 dB. For a JPA emitting a 
squeezed state and a 50 Vt termination emitting vacuum 
fluctuations into the input ports of our hybrid ring, the 
state in the beam splitter output ports is expected to be 
path-entang lecP. 



C. The cross correlation detector 

In Fig. a simplified sketch of the cross correlation 
detector is shown. Along each path, the signal emerging 
from the hybrid ring is linearly amplified, filtered, down- 
converted to an intermediate frequency /if = 11 MHz, 
and digitized at a sampling rate of 150 MHz by 16 bit 
analog-to-digital converters. The IQ-mixers used for 
down-conversion are biased by a strong local oscillator at 
5.626 GHz and split each output signal into its in-phase 
(Ii,2) and quadrature (Qip) components - therefore four 
ADCs are required. The digitized signals are finall y fed 
into an FPGA logic (details can be found in Sec. II A) 
which computes all correlations up to the fourth mo- 
ment in amplitude in real time. Phase synchronization 
is guaranteed by using a joint local oscillator for down- 
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Fig. S4. Sketch of the cross correlation detector. Colored 
arrows: output signals from the hybrid ring. Triangular sym- 
bols: amplifiers. Circles with crosses: IQ-mixers. 



D. Detailed setup 

A detailed setup of the experiment is shown in Fig. ^5| 
JPA, 50 Q terminated hybrid ring, and measurement cir- 
culator are anchored to the base temperature plate of a 
dilution refrigerator, whose temperature is stabilized to 
50 mK measured on the JPA sample box. Near the 50 Q 
termination of the hybrid ring, we measure a tempera- 
ture of 40 mK. The JPA signal and pump lines are heavily 
attenuated at various temperature stages. The coldest 
attenuator of the signal line is only weakly coupled to 
the base plate and the lower step exchanger of the fridge. 
Its temperature can therefore be controlled in the range 
^att = 50 — 800 mK, while all other components retain a 
stable temperature. This attenuator constitutes a broad- 
band microwave black body emitter, which is used to 
calibrate the gains of the amplification paths in a Planck 
spectroscopy experiment)^. The total power of each am- 
plification path detected at the ADCs is 



<i?,2> + <0?,2) 



R 



2 V2fcBT at t 



'1,2 



, (1) 



where R = 50Q is the input resistance of the ADCs 
and k B = 1.38 xlO -23 J/K the Boltzmann constant. The 
product of the gain G±^ and the photon number conver- 
sion factor K J = Rx2xBW mter xhf = 1.83xlO- 16 V 2 re- 
lates the measured autocorrelations (i 2 2 ) and {Q\ 2 )i 
which have units of V 2 , to the number of photons of fre- 
quency fo referred to the attenuator. The gain and 
the number of noise photons added by each amplification 
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Fig. S5. Detailed experimental setup. The bold numbers indicate the temperatures of the corresponding colored boxes. The 
two amplification and detection paths are marked with magenta color. The colored arrows denote the path entangled output 
state of the beam splitter. 
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path, ni 5 2, are retrieved from fitting Eq. ([!]) to the experi- 
mentally recorded autocorrelations. For pedagogical rea- 
sons, all formulas given in this manuscript assume equal 
gains and losses for the I and Q branches within each 
path. However, in the actual evaluation of the data, we 
do not make this assumption. We model the losses with 
beam splitters and also account for temperature gradi- 
ents along our cables. In this way, we make individual fits 
for the dependence of (Jj 2 ), (J|), (Qi), and (Ql) on T^tt- 
As an example, we show data and fit for (I±) in Fig. ^6] 
We first note that from this figure, we immediately see 
that the number of thermal photons in the mode fo is 
negligible at 40 — 50 mK. Furthermore, with a total loss 
of 1.8 dB between attenuator and beam splitter input, we 
retrieve Gdi/2 = 116.5 dB and n^i = 24.3 photons. Here, 
the index "d" denotes that Gdi/2 and n^i are referred 
to the input of the hybrid ring. From our reference-state 
analysis described in Sec. |IIIB[ we obtain, with respect 
to the beam splitter output ports, noise temperatures of 
3.00 K and 3.27 K for the two amplification paths. Con- 
sidering that our beam splitter reduces the input signal 
by 3.5 dB, the value of 3.00 K is in very good agreement 
with that of n^i quoted above. 




these elements in our discussion here since they serve 
purely technical purposes (e.g., debugging or switching 
to other experiments) and are not essential for the func- 
tionality of the cross correlation detector and the entan- 
glement detection. 



II. PROTOCOL FOR CROSS CORRELATION 
MEASUREMENTS 

A. Data acquisition 

In our experiments, we use either a squeezed state or a 
coherent state at one of the beam splitter inputs, while a 
vacuum state is always incident at the second one. Dur- 
ing measurements on squeezed states, the JPA pump at 
frequency 2fo is operated in pulsed mode. Similarly, we 
also use pulsed probe signals when measuring coherent 
states of frequency fo. The rise and fall times of the 
pulse envelopes are approximately 10 ns each. As shown 
in Fig. the measurement window always contains an 
off-region as a reference in addition to the signal. At low 
enough temperatures, this referen ce sta te can be consid- 
ered as the vacuum (see also Sec. IIIB). 




200 400 600 800 

T att ( mK ) 



40 60 
Time (|js) 



Fig. S7. Measurement protocol for squeezed and coherent 
states. During the high time ("1") of the pulse envelope, 
either the coherent signal or the JPA pump tone is on, during 
the low time ( "0" ) both are off. The shaded area denotes the 
measurement window. 



Fig. S6. Gain calibration. Dependence of the second moment 
(Ii) on the temperature of the attenuator temperature T a tt- 
Black symbols: data. Red line: fit. 

Apart from the key functional elements described in 
Sec. |I C[ the amplification paths contain isolating ele- 
ments at various temperature stages to avoid spurious 
correlated noise contributions. Furthermore, tunable at- 
tenuators after the IQ-mixers allow for a prebalancing of 
the four channels. The JPA pump and signal microwave 
sources, the local oscillator for the IQ mixers, and the 
clock for the ADCs and the FPGA logic are synchro- 
nized with a Rubidium-based 10 MHz source. The vector 
network analyzer is used to measure the JPA operating 
point, signal gain, and idler gain. Finally, in addition to 
the elements shown in Fig. £[5j several components such 
as mechanical microwave switches, power dividers, and a 
spectrum analyzer are used in the real setup. We omit 



For the dual-path reconstruction and the entanglement 
detection, the orthogonal / and Q quadratures of the 
noisy signal have to be recorded for both paths and prod- 
ucts of the type (l{l%Q™Q%), where j + /c + ra + n<4 
and j, k, m, nG No, need to be calculated for each 
recorded data point. However, the presence of the ampli- 
fier noise of our paths requires significant averaging. As 
a consequence, data transfer rate and computation time 
become a serious bottleneck in a computer-based acqui- 
sition system. By streaming the data from the ADCs 
directly into an FPGA logic, we solve these problems 
and are able to perform the moment calculations in real 
time. 

Data acquisition is triggered every 100 ills. For each 
trigger event, 8192 consecutive points are digitized in 
each of the four channels at a rate of 150 MHz. As shown 
in Fig. this results in a duty cycle of 54.6%. The 
data is streamed directly into the FPGA, where it first 
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passes a gain balancing and a digital down-conversion 
stage. The latter also allows for fine-tuning of the phase 
difference between the two paths to 180° and for correc- 
tion for phase imperfections of the IQ-mixer. In the next 
step, the number of data points is reduced to 512 with 
a digital cascaded integrator comb (CIC) filter. The fi- 
nal bandwidth is determined by the subsequent digital 
finite impulse response (FIR) filter. These filters also de- 
termine the measurement bandwidth (noise bandwidth) 
of 2 x BWmter = 2 x 489 kHz of the cross correlation de- 
tector. The corresponding transmission characteristics is 
shown in Fig. ^8j We note that one of the key advan- 
tages of digital filtering is that the frequency dispersion 
is flat and that the transmission characteristics are ex- 
actly equal for all four detection channels. Furthermore, 
for JPA signal gains up to 20 dB, the JPA bandwidth is 
larger than the measurement bandwidth. 
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Fig. S8. Measured transmission characteristics of the cross 
correlation detector. 
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Fig. S9. Typical averaged time traces of selected second 
moments from a squeezed state measurement with 0° phase 
and for 10 dB signal gain. Each averaged trace consists of 
5xl0 5 single traces. The rise time of 650 ns is determined by 
the digital filters. The step between JPA pump off (vacuum) 
and on (squeezed state) is shifted by 4 (xs with respect to the 
pulse shown in Fig. 371 because of a delay due to filtering. 



and amplitude a=(a) extracted from the reconstructed 
moments against power P gen at the output of the signal 
generator. Here, a) and a are the field operators of the 
input state as defined in Sec. Ill A| The expected lin- 

AP gen and 



ear and square root dependences, n(P gen ) = 
a (Pgen ) = ByJ P gen , are clearly reproduced. Within an 
error bar of less than two percent, B is the square root 
of A for independent fits. 



Finally, the moments up to fourth order in amplitude 
are calculated. For each moment and each data point, 
the average over a specified number of trigger events, in 
the following called ensemble average, is stored inside the 
FPGA. Figure ^9] shows typical time traces for selected 
second moments of a squeezed state averaged over 5x 10 5 
traces. Since the phase angle in this example is chosen to 
be 0°, (Qi) increases significantly above the vacuum level 
and (if) decreases below the vacuum level as expected. 
The cross moment (Iil 2 ) shows the characteristic cancel- 
lation of the uncorrelated noise contributions of the am- 
plification paths, which lies at the heart of the dual-path 
tomography. A Wigner function reconstruction based on 
this dataset is shown in Fig. fjnjb). 



B. Dual-path tomography of coherent states 

We first test the dual-path setup against coherent 
states of frequency fo = 5.637 GHz from a microwave sig- 
nal generator. The JPA pump is off in these measure- 
ments. In Fig. £ 10 , we show the photon number n= (a) a) 




P (mW) 

gen v ' 



Fig. S10. Coherent state reconstruction. Photon number n 
and amplitude \a\ as a function of the signal generator power 
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C. Phase stabilization protocol 

The phase stability between signal/pump and local os- 
cillator is better than ±0.3° for 5xl0 5 traces. Hence, re- 
constructions based on such a number of averages exhibit 
very good phase control as shown in Fig. £11 However, 



for quantifying the path entanglement properties, an av- 
erage over 8xl0 6 — 3xl0 7 traces is necessary to reduce 
the influence of the noise added by each amplification 
path to a negligible amount. In such measurements, the 
phase stability of our setup is not sufficient. For this 
reason, we record the data in 5xl0 5 trace portions and 
adjust the relative phase between signal/pump and local 
oscillator in a way that phase drifts are compensated. In 
particular, the data shown in Fig. 3 of the main text are 
recorded in this fashion. 
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Fig. Sll. Phase control for squeezed state reconstruction. 
Wigner functions of (a) the vacuum, (b) a squeezed state at 
0°, (c) a squeezed state at 1°, (d) a squeezed state at 45°, 
(e) a squeezed state at 90°, and (f) a squeezed state at 135°. 
The number of averaged traces is 5 x 10 5 . The residual ther- 
mal population of the vacuum is 0.102±0.005 photons. For 
the squeezed state, the JPA signal gain is 10 dB and the re- 
constructed photon number n = 8.67 varies by approximately 
0.5% for the different phase angles. 



III. THEORY: DUAL-PATH TOMOGRAPHY 
AND ENTANGLEMENT DETECTION 

In the optical domain, efficient single photon detec- 
tors and optical homodyning are established measure- 
ment techniques for quantum correlations^. However, 
despite recent effort J 37 " 39 *, the translation of these meth- 
ods to the microwave regime remains difficult because 
of the low photon energy. Therefore, quantum mi- 
crowave reconstruction based on off-the-shelf noisy lin- 
ear amplifiers and cross-correlation techniques, the cross- 
correlation detect oJ^, was proposed and experimentally 
implemente d 40 * 41 -*. Later, a second reconstruction tech- 
nique using only a single amplification path and decon- 
voluti on bas ed on a reference measurement was demon- 
strated 4 ^ 4 ^. In what follows, we describe the mathemat- 
ical details of these methods, which we a adopted to the 
needs of our experiments: the reconstruction of the state 
incident at the input port and the detection of entangle- 
ment between signals propagating in the output paths of 
the beam splitter. 



A. Dual-path reconstruction of the input state 

For the input state reconstruction, we can make use 
of the beam splitter and cross correlations. As described 
below, we in this way cancel the amplifier noise obscur- 
ing the signals because the noise contributions of the 
two paths are independent. The functionality of mi- 
crowave beam splitters is well-established for classical sig- 
nals a nd was r ecently demonstrated also for the quantum 

•vi 40141 1 

regim e 1 1 '. 

For quantum microwaves, we need to take into account 
the orthogonal signal quadratures I x 2 and Q x 2 . These 
are real valued voltages measured at the outputs of the 
IQ-mixers (see Fig. ^4] and Fig. ^5|. We can now define 
the dimensionless complex envelope functions 



fl,2 = ( J l,2 +*Ql,2)A/« J 



(2) 



where n is the photon number conversion factor intro- 
duced in Sec. |ID| The corresponding operators, 



(3) 



£l,2 = (A,2 +*Ql,2)/V« J 

can, in this situation, be expressed aPS 

£l,2 = ^1,2 + ^1,2 ' 



(4) 



Here, is the bosonic annihilation operator of the 

input signal into the IQ-mixer. The noise added by 
the latter is represented by its bosonic creation oper- 
ator v\ 2 . We now have [Ci 2 5 ^1 2] = 0? anc ^ the corre - 
lations (i{(ilr'il\il) n ') = '((il) m 'i((ii) n 'il') can be 
expressed in terms of the the measured (I1I2QTQ2) vldi 
Eq. hj by identifying £ 12 with £ 12 and £} j2 with £* ?2 . We 
note that, while in general, f, k' ,m' , ??/, j, /c, m, n £ No, in 
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this manuscript we restrict ourselves toj + k + m + n<4 
or, equivalent ly, f + k' + ra' + n' < 4. 

Using the beam splitter relations and the standard 
quantum model for linear amplifiers 45 , we can now write 



e 



i = y% ( + a + *) + v^d^iM + *I ( 5 ) 

(6) 



G 



( - a + v) + 



d2 ± ^2 



for our setup. Here, a and -0 are bosonic annihilation op- 
erators. They describe the modes incident on the signal 
and the 50 fl terminated input port of the beam splitter, 
respectively. The noise fields added by each amplification 
path are represented by the bosonic creation operators 
h\ 2 . Their effective temperature is mainly determined 
by the noise temperatures of the cold HEMT amplifiers 
and the cable losses between beam splitter and HEMT 
amplifier. Gdi,2 are the gains of the output paths cali- 

In order to simplify the 



brated as described in Sec. ID 



notation, we also define the operators 



$1,2 



G 



dl,2 

2~~ 



(y/ G di,2 -^1,2 + «i, 2 ) (?) 



u dl,2 



(8) 



We note that Vi^/V?- is a bosonic operator, as 
v 

["72"' "72"] = ^' ^ n ^ s wa ^' we arr ^ ve a ^ the simple ex- 
pressions 



Si 



■ a + v + V{ 
■a + v + Vj . 



(9) 
(10) 



We note that the operators Vi 5 2, a, and and therefore 
also 5i 5 2, are referred to the input of the beam splitter. 
With these definitions, we can generalize the dual-path 
state reconstruction technique, which we originally de- 
veloped in Ref. |40j in a way that it becomes applicable 
to the experimental setup presented in this work. Mak- 
ing the reasonable assumptions that v is a weak thermal 
state with a measured temperature of 40 mK and that 



(V 1 ) = (V 2 ) = 



(11) 



for the noise added by the amplification paths, we recur- 
sively obtain the signal moments 



l t l — l-L m 1 m—m 1 — l l 1 —k 1 l — l 1 —k 2 m \—3\ m — m 1 — j 2 x /7 j \ / \ / 
fc 1= /c 2 =0.^=0 j 2 =0 fcj_=0 fc 2 =0 j(=0 j 2 =0 

^i - M - k 2 \ ( m \ - 3i\ fm-m 1 -j 2 



k[ 



Ji 



J2 



^1/ V ^2 J \Jl J \ J2 

^_^l-l 1 +rn-rn 1 +j 2 +k 2 



x ((fi + ) fc i +fc a a^+^)(({)t) fe i+ /c 2 gii+i^^i- fc i- fc i (T> 1 t ) m i-^'i-^)(T> 2 

i x i— m 1 — l l^—ki l—l 1 — k 2 m\—3\ 

EE E E E E 



y-Z Z x k 2 k 2 /VVf 



^1/ V ^2 J \3i J \ 



/bo 



fc 1= ofc 2 =o j 1= o fe^=o fc^=o j;=o 

x (-ly-'i+fca^fit^+fca a^ +m - m i)((t) t )' s i + ^ ^)(t^ 1_fel ~* i (y i t) TO i-^-^)(y 2 '- i i-^-^ 

E ' E 1 (I 1 ) f ^0 (-i)'-' 1+fe2 ((^^ a™>(^>(^) 



) 



fc 1= A; 2 =0 



^ =0 



E (MW 1 *-' 1 *")^ 1- * 1 ) 



(12) 



for Z,m,Z 1 ,m 1 G No from the measured noisy correlations. In this process, we also have to compute the noise 
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moments 



(v{(vir) = {(sirst) 

r s — 1 r—ki s—Ji 

-EEEE 

fc i=° 3i=° k 'i=° 3i=° 
r-1 



r — k x 



3 )<(a t )* 1 o^X^t)--*!-*! ^-^-^Xy^^)^) 



(13) 



<w 2 t> = <(s 2 ysi> 
-EEEE 

fc i=° ii =0 fc i=° ii=o 

r-1 



r \ / s\fr — k x 
kj\jj\ K 



E ( fc r i )(-i) r " fei ((« t r fei )(^ i (^ t ) s > 



5 /^(-^^^(^^^^^(^^""^"^^"^"''^(^u^F 1 ) 



(14) 



/c 1= 

associated with both amplification paths for r + 8 > 1 and r, s E No, again in a recursive fashion. 



In other words, the formulas for the moments of order 
I + m are established using those of the moments of order 
Z + 771 — 1. The formulas obtained in this way are not 
unique, they depend on the specific choices of l x and m 1 . 
We find that the statistical uncertainty in our results is 
minimized by using the mean value of all formulas found 
for constant I + m. 



From the reconstructed signal moments of Eq. (12), 



we can readily extract the squeezing below the vacuum 
in decibel, 

101og 10 ( - (a 2 )e-^ - ((a + ) 2 )e^ + 2(a) a) + 1 

+ (a) 2 e-^ + (a^V* - 2(a t )(a)) . (15) 



I 

Here, 



the 

<d> 2 = 



angle 

I (a 2 )- 



(j) is defined via the relation 
(a) 2 |e 2 ^ and the argument of the log- 



arithm is the ratio between the variance of the squeezed 
quadrature and the vacuum variance. 



In the case of infinitely many reconstructed moments 
((a^) l a m ), the Wigner function W(q,p) of an arbitrary 
state can be completely reconstructed. However, in this 
work we record these moments only up to fourth order, 
I + m < 4. As we find that the higher moments are con- 
sistent with those of Gaussian states, we are allowed to 
restrict ourselves to moments with I + m < 2. This en- 
ables us to use an analytical approach ! 46 * 47 !, which yields 



W(q,p) 



: exp 



{y + 1/2)|C - (d)| 2 - (M*/2)(C ~ (a)? ~ WgKC ~ (« f » 2 

1/2) 2 -M 2 



(16) 



with ( = q + ip, /i= ( a 2 ) - (a) 2 , and v = (a) a) - \(a)\ 2 . 
As explained in Sec. |III A[ we have chosen our defini- 
tions such that phase space variables q and p are di- 
mensionless and their value represents the square root 
of a photon number. Since any Gaussian state can be 
written as a displaced squeezed thermal state, we can 
also extract the effective mode temperature anal}/ 
from the reconstructed input state moment s^ 6 -^ 
note that this effective temperature contains contribu- 
tions from the physical temperature and from losses. 



ilytically 

m we 



B. Reference-state analysis of the output state 
moments 

In order to detect the entanglement between the two 
paths independently from the dual-path reconstruction of 
the input state, we cannot assume that the hybrid ring 
is working as a beam splitter. We therefore follow a dif- 
ferent route: we reconstruct the moments of the output 
state by means of a calibration against a well-known ref- 
erence signaP^EU. The obvious choice for this reference 
signal is the two-mode vacuum. In this way, the beam 
splitter is treated as a black box device which, for a vac- 
uum state at each input, produces uncorrelated vacuum 
states at each output port. This rather general assump- 
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tion holds well for the temperatures measured for atten- 
uator and termination, 40 — 50 mK. In this situation, the 
complex envelope operator becomes 



G rl 8 1 + y/G rl 
\A*r2 




lft 



lft 



(17) 
(18) 



Here, s x 2 is referred to the output of the beam splitter, 
and G r i,2 is the effective gain of the amplification paths. 
Note that the G r i,2 are numerically different from Gdi,2 
because they do not contain the beam splitter losses. Af- 
ter defining the operators 



V lt2 
^1,2 



2 1 ^1,2 + ^1,2 



~£l,2 ' 



we again arrive at the simplified expressions 
S 1 = s 1 + V? 



S 2 — s 2 



Vr 



2 ' 



(19) 
(20) 



(21) 
(22) 



We now evaluate the correlations of the outputs of the 
two channels, 

((stM mi (4 f )' 2 4 m2 > 

= <(4 + ^i)' 1 ^ + vir^si + v 2 y*(s 2 + v}r*) 

l ± l 2 m-L m 2 

= EEEE 

fe 1= fc 2 =0 j 1= j 2 =0 



Ji 



32 



x ((at)'i-*i^i-Ji(4)'2-* 2 ^-J a ) 



(23) 



With the terms ((4)' 1 fel Si 



2 s 2 2 2 ), which 



can be calculated straightforwardly for our reference 



forms 



a system of linear equations, 
ows us to extract the noise terms 

Uj 2 \ 



state, Eq. 
The latter al 

(V?KVi) jl V£ 2 (V2) j *) related to our amplification paths 
by algebraic inversion. Once knowing these noise terms, 
we can extract the signal correlations for the squeezed 
state input again from Eq. (23) and algebraic inversion. 
We note that also more sophisticated reference states and 
device models, such as thermal states incident at a beam 
splitter, can be chosen as reference state. Depending 
on the pre-characterization of the used components, this 
approach might account better for experimental imper- 
fections. Although the latter typically tend to reduce 
the degree of entanglement, our entanglement detection 
turns out to be quite robust against them. 

Altogether, the method described above would, in 
principle, allow for a reconstruction of the output state as 
far as this is possible with four moments. We note that 
from such a reconstruction also all entanglement prop- 
erties could be derived. However, as shown in the next 
section we choose a different approach. 



C. Entanglement witness 



The detection of entanglement requires substantially 
less information than a full state reconstruction. Indeed, 
there exists an infinite number of witnesses and criteria 
which allow one to decide whether or not a state is en- 
tangled. Each of these criteria uses only a small amount 
of information about the examined state. In this work, 
we use the witness matrix 



M« EE 



/ 1 (sj (s\) 

(4) (4*i) ((4) 2 ) 

(h) (si) i + (s\h) 

V(4) <*i4> (44) 



(4) 
(44) 

(h4) 



(* 2 ) 

<4«2> 



\"2l 



(*2*2> <*2> 
{(si) 2 ) l+(s\s 2 )) 



(24) 

which contains up to second order moments of the beam 
splitter output state. If M^ 2 ^ has at least one negative 
eigenvalue, the state is entangled. The absence of a neg- 
ative eigenvalue implies separability only in the case of 
Gaussian states. 



D. Negativity 



For a bipartite system, the amount of entanglement 
between the subsystems A and B can be quantified by 
means of the negativity 



(25) 



where p is the density matrix of the total system, and 
\\p TB ||i = Tr\p TB | is the trace norm of the partial trans- 
pose of p with respect to subsystem B, p TB . If M(p) > 0, 
the state is entangled. For a maximally entangled state, 
M{p) ->>oo. 

In the case of Gaussian states, all measures of entan- 
glement are equivalent, and they are defined by the co- 
variance matrix 



ol 7 

7 T 



(26) 



Here, we define the matrices 



«i as) g= ( Pi & 



In 7i2 

721 722 , 



(27) 
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with 

Oil 

a 2 
a 3 

Pi 
ft 

7n 



721 



722 



(s 2 ) + ((s\) 2 )+2(s\s 1 )-(s 1 +s\) 2 + l 
-(s 2 1 )-((s\) 2 )+2(s\s 1 ) + (s 1 -s\) 2 + l 

i(-(si)+((s\) 2 )+(h) 2 -(m 

(si) + ((4) 2 ) + 2(4s 2 )-(s 2 + 4) 2 + i 
-<*!> - ((4) 2 ) + 2(4^) + & - 4) 2 + 1 
K-(4) + ((4) 2 ) + (h) 2 -(4) 2 ) 



+ (S 2 S! + S 2( sJ + 4«1 + 44)/ 2 

- (si + 4)(^ + 4) 



712 = (SiS 2 - S 1 S 2 + S l 5 2 - ^I4)/ 2 ^ 

+ (^Si + S 2 s| - 4«1 - 44)/ 2i 

+ 4)(^2 -4) 



(s 1 S 2 + ^1^2 4^2 44)/ 2 * 



+ (Mi 



^2^1 



n D 2 

■ 4-si 



at at 



Sl>/2i 



+ -4)(^2 + 4) 

+ (-Mi + Ml + 4^i - 44)/ 2 

+ (h- 4>(^2-4> 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



(34) 



(35) 



(36) 



(37) 



and 7 T being the transpose of 7. Finally, the negativity 

becomes^! 



E. Verifying consistency with a Gaussian state 
using higher order cumulants 

In order to check whether the states we reconstruct are 
consistent with Gaussian states, we evaluate the (I + m) th 
order cumulants (((a^a m )) for / + m<4 and Z,raE No- 
Equivalent ly to the moments, the cumulants describe a 
probability distribution. The definition of cumulant for a 
quantum state with density matrix p can be written 

a*P 



gl gm 

d(ip*y d(ipy 



2 



lnTr(e^* at+ ^ a V) 



J/3,/3* 



=0 

(39) 



Only Gaussian states have a finite number of nonzero cu- 
mulants. More specifically, all their cumulants vanish for 
I + m > 2. In other words, finding a nonzero cumulant of 
3rd or higher order implies that the state is not Gaus- 
sian. Despite not being a strict proof, the fact that the 
3 rd and 4 th order cumulant are very small or vanish in 
an experimental reconstruction constitutes a reasonable 
indication that the reconstructed state is Gaussian. 

As an example we spell out the 3 rd order cumulants 
as functions of the moments at the beam splitter outputs 
reconstructed with the reference-state method. We find 52 



<<sl2>> H^2> - 3<<2><Sl,2> + 2<S 1>2 }^ 



(40) 



«4,2^2» = (4 



2*1,2/ 



<4,2><*1,2> 



2 (4,2^ 2 )(^2) + 2(4o)( S V 2 ) 2 . (41) 



:|o,^J =max{0,A/-}, (38) 



where v = y (a(<t) - y/A 2 (<r) -4 deter) /2 and 
A(cr) = det a + det /3 - 2 det 7. 

Note that, despite not being a measure, the negativity 
kernel J\f is a witness for arbitrary bipartite entangle- 
ment. In fact, if a non-Gaussian state has the same 
first and second moments as an entangled Gaussian 
state, it is entangled 50 . Consequently, J\f>0 implies 
entanglement for any bipartite state. 



ACKNOWLEDGMENTS 

The authors thank Christopher Eichler for discussions. 
We acknowledge support from the Deutsche Forschungs- 
gemeinschaft via the Sonderforschungsbereich 631, the 
German excellence initiative via the 'Nanosy stems Initia- 
tive Munich' (NIM), from the EU projects SOLID, CC- 
QED and PROMISCE, from MEXT Kakenhi "Quantum 
Cybernetics", the JSPS through its FIRST Program, the 
Project for Developing Innovation Systems of MEXT, the 
NICT Commissioned Research, EPSRC EP/H050434/1, 
Basque Government IT472-10, and Spanish MICINN 
FIS2009-12773-C02-01. 



* Edwin . Menzel@ wmi . bad w-muenchen . de 

(1 T. Yamamoto, K. Inomata, M. Watanabe, K. Matsuba, 

T. Miyazaki, W. D. Oliver, Y. N akamura, and J. S. Tsai, 

|Appl. Phys. Lett. 93, 042510 (200871 
' 2 R. E. Collin, Foundations for Microwave Engineering, 2nd 

ed. (Wiley and Sons, New York, 2001). 
13 M. Mariantoni, E. P. Menzel, F. Deppe, M. A. Araque Ca- 

ballero, A. Baust, T. Niemczyk, E. Hoffmann, E. Solano, 



A. Ma rx, and R. Gross, |Phys. Rev. Lett. 105~ 33601 
(2010)| 

34 E. Hoffmann, F. Deppe, T. Niemczyk, T. Wirth, E. P. 
Menzel, G. Wild, H. Huebl, M. Mariantoni, T. Weifil, 
A. Lukashenko, A. P. Zhuravel, A. V. Ustinov, A. Marx, 
and R. Gross, [Appl. Phys. Lett~97 , 222508 (2010)] 

35 Z. Y. Ou, S. F. Pereira, H. J. K imble, and K. C. Peng, 
Phys. Rev. Lett. 68, 3663 (1992)] 



16 



36 U. Leonhardt, Measuring the Quantum State of Light 
(Cambridge University Press, Cambridge, 1997). 

37 B. Peropadre, G. Romero, G. Johansson, C. M. Wilson, 
E. Sola no, an d J. J. Garcfa-Ripoll, Phys. Rev. A 84, 



063834 (2011) 

^ Y.-F. Chen, D. Hover, S. Sendelbach, L. Maurer, S. T. 
Merkel , E. J. Pritchett, F. K. Wilhelm, and R. McDer- 
mott, |Phys. Rev. Lett. 107, 217401 (2011)] 

39 A. Poudel, R. McDerm ott, and M. G. Vavilov, 
|arXiv:1201.2990vl (2012)] 

40 E. P. Menzel, F. Deppe, M. Mariantoni, M. A. Araque 
Caballero, A. Baust, T. Niemczyk, E. Hoffmann, A. Marx, 
E. Sola no, and R. Gross, [Phys. Rev. Lett. 105, 100401| 
(2010) 



41 D. Bozyigit, C. Lang, L. Steffen, J. M. Fink, C. Eich- 
ler, M. Baur, R. Bianchetti, P. J. Lee k, S. Filipp, M. P. 
da Silva, A. Blais, and A. Wallraff, |Nature Phys. 7, 154| 

I (20H)| 



42 C. Eichler, D. Bozyigit, C. Lang, L. Steffen, J. Fink, and 
A. Wallraff, [Phys. Rev. Lett~1 06, 22050 3 (2011)| 

43 C. Eichler, D. Bozyigit, C. Lang, M Baur, L. Steffen, J. M 
Fink, S. Filipp , and A. Wallraff, |Phys. Rev. Lett. 107,| 



and A. Blais, 



113601 (2011) 

M. P. da~Sn7a, D. Bozyigit, A Wallraff 
|Phys. Rev. A 82, 043804 (201 0)] 

45 C. M. Caves, |Phys. Rev. D 26, 1817 (1982)1 

46 V. Buzek, G. Adam, and G. Drobny, |Phys. Rev. A 54, 
804 (1996) 

V. Buz ekTG. Adam, and G. Drobny, |Ann. Phys. 245, 37 
(1996)| 



E. Sh chukin and W. Vogel, Phys . Rev. Lett. 95, 230502 
(2005) 



G. Adesso and F. Illuminati, Phys. Rev. A 72, 03 2334 
(2005) 



P. Hyllus and J. EisertjNew J. Phys. 8, 51 (2006) 



R. Schack and A. SchenzleT lPhys. Rev. A 41, 3847 (1990) 
L.-Z. Jiang, Proc. of SPIE 7846, 784612 (2010) - 



